Constraints in the coupling Star-Life 



Luigi Sertorio ^ 



Giovanna Tinetti ^ 



Dipartimento di Fisica Teorica - Universita di Torino. 
Istituto Nazionale di Fisica Nucleare - Sezione di Torino 

DFTT 27/2001 



sertorioOto. infn.it 



tinetti@to.infn.it 



Abstract 



If life is sustained by a process of photosynthesis, not necessarily the same existing on Earth, 
the surface temperature of the star and the orbit of the host planet cannot be whatsoever. 
In fact the global life cycle, no matter how complicated, must contain in general an upper 
photochemical branch and a lower dark branch, characterized by a higher and a lower 
temperature (§2). These two temperatures are star-orbit related. The velocity along the 
cycle or, in other words, the power of the life machine, depends in general on several other 
parameters. First of all the Gibbs photon availability, which is a star-orbit parameter and 
is the input for the upper branch (§3). Then follows the energy cascade that develops along 
the organic web with a large number of interactions and typical times that must match the 
typical times generated by the combination of spin value and orientation, eccentricity and 
precession (§4). Finally, the capacity of the web to keep the global life cycle running along 
the life span of the star, comes from some inner form of self-endurance and self-balance. 
The property of not being transient could be the right way of introducing the concept of 
intelligent life (§1). 
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1 Life from inside 



We need to face a fundamental distinction, namely among deterministic and finalistic 
systems. 

Deterministic systems are described by equations of motion, either known ex- 
plicitly or supposed to exist, which come from a minimum principle related to a 
Lagrangian. This is the domain of physics. 

Finalistic systems are those in which what is minimized is another functional 
that can never be reduced to a Lagrangian type, because the object to be minimized 
contains "external inputs". This is a formal way of introducing the concept of 
finalism and can be helpful, although not necessarily the profoundest. 

A finite size finalistic system is realized by a standing chemical structure, occu- 
pying a certain domain in space. Such standing finite structure requires a power 
density. This object can be called organism, including with this word a large variety 
of structural complexities and control abilities. Going on in the process of consider- 
ing complex structures and related complex control functions, we cross a borderline 
and talk of living organisms. 

Finally we wish to touch only marginally the concept of intelligent life. This 
concept is particularly ambiguous, because the only way we know to communicate 
acquired knowledge is a forward process both for an individual and for a global com- 
munity. Even the most abstract way of thinking, mathematics, develops in time. 
Therefore the question is: intelligent at what time, rather than intelligent up to a 
certain score. It follows that for a given realization of life, for instance that on Earth, 
or another, rather than abstract intelligence, we should talk of history of intelligence 
assessed in a particular time interval of the presumed course of development. This 
also implies typical times. For instance the human terrestrial civilization is under- 
standable for, say, few thousand years. This may be a very long time interval for 
another form of life on a different planet, or extremely short, almost not noticeable, 
for another. In all cases, it is the global cycle, and its history, that contains the 
maximum amount of intelligence that was given to that particular planet. 

In conclusion, the discussion of intelligent life referred to another form of life 
existing elsewhere, remains a rather uncertain issue. 
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2 Life from outside 



The analytic knowledge of a living organism assumes, as a methodology, the con- 
sideration of a simple component, or elementary organism, then several components 
constituting a higher organism, and so on. Finally animals, and man, supposed to 
be the best, or highest organism. Let us go on, a step further, and consider the 
whole of the living system, however complicated. This is what is seen from outside. 

The complex operation of this "global organism" is peculiar. In fact the power 
input is given uniquely by the incoming photons and the power output is again 
only by photons, those leaving the planet. This is the way in which global life is 
organized, if it is photon life. So while every living organism has an ambient, an 
ecosystem around it, the global organism is that one having for ecosystem the entire 
Universe, appearing via the star nearby and the cosmological blackbody background 
radiation. 

We have considered in a preceding paper IQ the case of a generic life based on 
a process of mining energy and material, and based on a process of extracting the 
thermal availability embedded in a thermal disequilibrium flow; this is the case of a 
planet without star, which slowly cools off and therefore contains thermal gradients. 
We belong to photon induced life. We wish to consider this case in general and not 
uniquely in the terrestrial realization. 

Photon life is a system open to energy in the form of in and out radiation fluxes. 
Concerning the material transit across the living web, we have the important con- 
cept of asymptotic constituents. Within the global organism a large number of 
complicated assemblies of molecules, organisms, are formed, stand for a while, and 
decay. In doing so, they communicate material and energy to each other. This is 
the global metabolism. Nevertheless, this global metabolism rests on the circula- 
tion of asymptotic molecules; those, among the released, which are capable of being 
photo-synthesized again. For instance, in the particular case of the terrestrial life, 
the asymptotic molecules are H2O and CO2. There is no global mining and corre- 
spondingly accumulation of products unable to return to the mine. 

The asymptotic molecules come out of the web at a certain rate, expressed in 
mass per second, move along the planet in either liquid or gas phase, and then return 
again as components in the process of photosynthesis to restart the construction of 
complex structures, which keep feeding the web, and so on forever, in principle. The 
global rate of photosynthesis must be equal to the global rate of dark release if 
the global circulation is stationary. It is in this sense that the molecules H2O and 
CO2 are not only participating in the dynamics, as any other component of the living 
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web, where nothing is fixed, but in addition are asymptotic. If the global life were 
terminated, all the molecules participating to the standing global structure would be 
released on the planet and would settle down to the final distribution, namely that 
corresponding to the minimum of the Gibbs potential, at those values of T, and P, 
which are pertinent to the planet. 

From the viewpoint of energy and entropy, the global web corresponds to this 
cycle in the diagram temperature-entropy (fig. |2.1| ). 



T 

T, 



Photosynthesis of H^O + CO, -> Polysaccarides 









Decay into HjO and CO, 



Figure 2.1: 



Ti is the Planck spectrum temperature of the photons activating the chlorophyll, T2 
is the dark temperature of the final decay which ultimately must reproduce H2O and 
CO2 in the right amount that is required to restart the cycle on the upper branch 
at Ti. 

We show in fig. the window of frequencies able to operate the photosynthesis 
in the case of the terrestrial carbon life. For another triplet H, C, O we may expect 
another Tg and another wavelength window. 

The energy balance implies Ei incoming with the photon flux originated by the 
star, E2 outgoing with the photon flux leaving the planet, which acts as the thermal 
bath at T2. The global cycle is stationary if 

El = E2 = E 

For the entropy S we have 
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Figure 2.2: 
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Let us consider now the concept of availability. The Gibbs availability [Q] 
spending to the standing structure which performs the cycle is 



corre- 




) 



Notice that we expect T2 <C Ti in general. At T2 the asymptotic molecules perform 
the decay reactions of the lower branch, or the dark branch. To fix the ideas, for the 
chemical life on Earth based on C, O, H, N, we have Ti = Tq ^ 5800 K, T2 ^ 300 K 
and 



Now, much more interesting than the above equilibrium thermodynamic description 
of the cycle, would be the understanding of the non-equilibrium thermodynamic 
description, that is the time rate of the various processes (see §3). From such an 
analysis, we could formulate the concept of power of the global ecosystem cycle. 
In fact the photons come as a flux, energy per unit time, per unit surface. They 
trigger all the dynamical interactions in the living web. These interactions have a 
large number of typical dynamical times. Consider for example the time pertinent 
to an electronic molecular quantum level excitations, r ~ 10~^^ s; and consider for 
instance, on the other extreme, the life time of a red tree 



AW ~ E 



r~ 2000years~6-10^"s 
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Between these two extremes, so many other typical times are excited and all of them 
must function combining with each other to produce the actual power balance 

bio bio 

where is the portion of the radiation flux coming from the Sun absorbed by the 
living ecosystem, and is the contribution of the ecosystem to the radiation flux 
emitted by the Earth. 

Another form of life may function with asymptotic components H, C, O, which 
are synthesized into another structure, a polymer 1^,(7,01, requiring a spectral 
temperature Tg and a dark temperature T. Different sets of typical times will be 
involved. But not whatsoever. We know that the laws of quantum mechanics must 
be obeyed, therefore we may think of a different chemical life but within the identical 
chemistry. 

The intermediary between a form of life and its star is the "protagonist planet" 
and its orbit. The typical times generated by the orbital parameters on one side, 
and the typical times of the unknown form of life, on the other side, must combine 
together. This is what we mean by "constraint" in the present paper. 
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3 The two temperatures and the year 



Consider now a planet that contains the appropriate set of molecules which can or- 
ganize themselves in complex chemical structures that eventually evolve into life. 
Among these molecules there must be the subset of asymptotic molecules, not nec- 
essarily H2O and C02- The minimum conditions for the existence of life are given 
by the existence of the photon temperature Tg and the dark temperature T. 

Notice that the life cycle dark branch implies not only a dark T but also a 
dark pressure P. The pressure P is a parameter that comes from the history of 
the formation of the planet itself. This history determines if there is an abundant 
atmosphere or not, and consequently the value of P. What is important is the 
constraint that the asymptotic molecules be in fluid state to ensure the necessary 
supply-flow of matter required by the incoming photons to perform their synthesis 
job. In fact the photosynthesis is a high energy density process and is fast; the 
corresponding dark decay is a low energy density process and is slow. To close the 
cycle a large surface is necessary in order to sustain the transport of matter towards 
the synthesis centers, hence the fluid motions. Roughly speaking we can think of 
this balance: 

^Tio = <yT^s, (3.1) 

where Si is the total area of the photosynthesis centers. Si can be thought, in the 
language of elementary particle physics, as the total cross section of the reaction 

7 + center — > inner operation on H2O, CO2 polisaccharide 

At the lower branch the organism decays into the asymptotic molecules that ther- 
malize at T and uniquely release the corresponding heat of the exothermic reaction 
into radiation with density of power a T^. The total out going flux is 

^"^^ = aT'S, (3.2) 

What is 53? Consider fig. 

Si is the area of the photosynthetic center 
^2 is the area occupied by living organism 

53 is the effective area needed by the asymptotic molecules to release their heat of 
reaction at temperature T 

An R'^ is the total emitting surface of the planet. The closure of the cycle implies 
$i" = 

^ mo ^ bio 
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Figure 3.1: 



namely 



aT^Si = aT^S3 



and this equation implies 



From appendix B, we have 



therefore /3 is dictated by astronomy. We expect 

to be a rather small number. For instance in the particular case of the Earth we 
have Ts = 6000 K, T = 300K and therefore 

/3 ~ 0.6 • 10"^ (3.4) 

An interesting question is: how large is ^3 with respect to 47r i?^? Reference |^] gives 
for 

$^^ = 10^5W 
so that we can write 

f7T^53 = 10^^W (3.5) 
On the other hand we know that 

4 vr i?2 cr T'' = 1.2 • lO^'^ W (3.6) 



Comparing (^) and ( |3.6D we get 

53 ~ 10-2 . 47ri?2 (3 7) 

The evaluation ( |3.4D indicates quite clearly that we need fluid transport of the 
asymptotic molecules; hence it makes sense to talk not only of a temperature T, but 
of an appropriate pressure P. 

How can we predict the temperature distribution over the protagonist planet? It 
is almost impossible to describe in general the fluid motions existing over an abstract 
planet. The best knowledge we have about the dissipative structures, comes from 
the simple and well studied case of the Benard cells. What we know qualitatively, 
in general, is that the presence of the fluid circulation increases the heat transport 
with respect to the transport described by the pure rigid planet, that is the heat 
transport described by the Fourier equation 0]. 

In conclusion we study a model planet which is a rigid black sphere. There are 
two motivations for this study. 

- The thermodynamics of a rigid surface is governed by the Fourier equation, which 

is linear, and this property allows a rather satisfactory analytic study of the 
temperature field. 

- The surface temperature distribution given by the rigid model is the most reason- 

able input for the actual, unknown, atmospheric dynamics and corresponding 
heat circulation. 

We can say that the black planet plays the role of the lower, warm, plate of Benard. 

Note: In the particular case of the Earth, the rigid black planet approximation 
is not used as an input for the understanding of the atmospheric circulation. This 
is because the Earth has a known geography, and we are sitting here to observe it; 
this knowledge is a much more detailed input for the study of the fluid circulation, of 
course. But this is the approach from inside and now. Different geography, different 
circulation. It is the historical anthropomorphic viewpoint. And this viewpoint is 
inadequate for our purposes. For the abstract planet which is the protagonist of the 
present study the black sphere model is the only model we can think of. 

We dedicate §4 to comments on the atmosphere. 

After these preliminaries, we consider a star with surface temperature Tg and a 
planet rotating around it in a closed elliptic orbit. Rigorously speaking, a two body 
planetary orbit is an idealization never realized in nature on the long time scales. 
Nevertheless this assumption is very attractive, because it allows us to push forward 
some general analytical results, so we adopt it. 
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The planet at position ■0 on the ellipse receives a total radiation flux (see ap- 
pendixes A, B) 



ttR^ 



(3.8) 



The photons come with Planck temperature T^. The function cP{ip) mod 2 7r is 
known. Let us assume that the incoming and outgoing fluxes are equal for any ijj 



aTf^ 



ttR^ = o-r4(i/;)4 7ri?2 



Notice this is a power balance not a linear relationship between temperatures. From 
this equation we can define a "radiation average" 



<r(^)>= T^{^) 



1 



Rs 



Rs 



V2^dW ' a(l-e2) 
The amplitude of the oscillation of < T{ip) > during the year is 

A < r >=< T > {1'^ solstice)- <T> (2""^ solstice) 

With our choice of the polar angle ip we have: 
ip{V^ solstice) = —5 and ^ip {2'^'^ solstice) = f , so we get from 



V^l - e sinV' (3.9) 



A < T >= 




1 



1 



^f2 V a 

Clearly this amplitude is zero for e = 0. 
The annual average of < T{;\p) > is 



1 



27r 
1 



2 7r 



< T > = — / < T{i,) > d^ = ^ 



Rs 



^2 V 2^ 



e sin ip dip 



R. 



V2 7r V « (1 + e) 



-2e 



4' 1 



^,37r -2e 



(3.10) 



E{(j),m) are elliptic integrals of the second kind known in the literature S. The 



dependence on e of < T > is obtained using the tabulation of E{(p,m), (fig. |3^ 
Finally we consider the ratio 



1 



< T >(e) _ l_ 

^T^(O) " vr V (1 + e) 



^,37r -2e 



4 ' 1 



1 



This formula comes from the observation that E{(p^ 0) = (p. 

Now, the life processes on the planet depend on the time t not on the angle 
ip (appendix A, eq. ([A.6|)). Our ellipse is characterized by a and e. We show in 
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appendix A that, given a and the masses of the star Mg and the planet nip, it 
fohows the annual period r^: 



Ty = 2TTj—— (3.11) 

At this point we observe that our planet cannot be of the Jovian type (mostly 
hydrogen). Consequently we neglect rup in ( 3.11 ) so that Ty depends only on a 



and Ms- But we know that for stars in the Main Sequence the mass Mg and the 
temperature Tg are related. Consequently if the hypothetical photosynthesis exists 
over the protagonist planet, that particular form of life implies a certain Tg and a 
certain Ti, (b stands for best approximation of the temperature of the dark branch). 
The couple {Tg, Tb) may couple only with one couple of {Mg, a); therefore the lenght 
of the year follows from the knowledge of {Tg, Tf,) 

'^y = Ty{Ts, Tb) 

Let us assume for simplicity that e = 0, d = a. We use for Tb the value ( |3.9| ), so Tb 
is a function of Tg and a, therefore given the couple Tg and Tb using ( p.llj ) we get 
Ty. We show in the table below Ty as a function of Tg and Tb, Ty is expressed in units 
of Tye, namely the terrestrial year. On the coordinate axis Tg we add the column tg 
(life of the star) and Mg (mass of the star) taken from B 
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0.5 
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1.67 


0.7 
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1.4 
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0.49 
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0.57 
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0.075 


9550 


18.18 


5.39 


2.27 


1.16 
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0.025 


12590 
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17.5 


7.36 


3.77 
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0.012 


15135 


120.5 


35.7 
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5-10-^ 
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310-^ 
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480.6 
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2-10"^ 


25120 


722.6 


214 


90 


46.2 








200 


300 


400 


500 



Tt{K) 



In the particular case of our carbon life (H, C, O), as it comes on Earth: 

Ts = Tq = 5800 K 
ts = tQ = 10^^ years 
e = 0.017 

R^ = Rq = 6.88 • 10^ m 
Ms = Mq = 1.99 • 10^0 kg 
a = 1.49 • 10" m ~ d 
nip = rrif. = 6 • lO^'* kg 
Tb = Tav.obs. ~ 288 K 
Tye = 3.15 • lO'^S 

From ( |3.9D we get 
< T >~ 279K 

The difference between the calculated < T > and Tav.obs. observed is explained 
by the combined effect of albedo and greenhouse due to the atmosphere. This ex- 
plains why in the column (300 K) and in the row Ms = Mq, we find a year a little 
shorter than the terrestrial year. 
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4 Atmosphere, yes and no 



Consider a protagonist planet equipped with the right amount of atoms, namely those 
that enter into the structure of the living organisms. Next consider the molecules 
that are the building blocks and which undergo chemical reactions all the time. In 
fact a living structure is different from the ensemble of the inert constituents. The 
inert, or lifeless, constituents stay in the physical-chemical state of minimum Gibbs 
potential pertinent to a given T and P, while the organic structure does not. Among 
the constituents we have the asymptotic molecules. For carbon life they are H2O 
and CO2. 

The photosynthesis cycle requires a small area occupied by the photosynthetic 
centers and a large area occupied by the asymptotic components (indicated by 5i 
and (S3 in §3). Therefore a process of release and return to the centers must be 
set into operation. This is called transport. Transport may come by diffusion and 
convection. Diffusion in solids is negligible; so we need asymptotic molecules in 
liquid and gas form. For carbon life on Earth CO2 is always gaseous, H2O sits in 
the three phases. Question: it is important, or even necessary, that at least one of 
the asymptotic molecules should explore all the phases? The answer is most likely 
yes. And the reason is the following. Both H2O and CO2 are greenhouse gases. 
Greenhouse gases are necessary to control both the average temperature and its 
oscillations. This will be shown numerically in §5. H2O creates greenhouse by direct 
molecular interaction with the photons, and creates albedo by indirect interaction. 
Greenhouse belongs to quantum dynamics, albedo belongs to classical optics. These 
interactions are linked by a triple feedback scheme as follows. 



In fact H2O evaporates from the liquid phase in bulk when the temperature T{f,t) 
and the pressure P{r,t) command it. So H2O moves to the atmosphere. Up there 
operates as a gas, and does the job of greenhouse. But the temperature and pressure 
fields in the atmosphere change from place to place and from time to time, so that 
the H2O sitting in the atmosphere may find itself in a local condition very close 
to condensation; it begins to develop a droplet, evaporates, condenses again and so 
on. That is the H2O in the appearance of cloud. A little cooler and the droplet 
becomes a drop and rains down; a little warmer and H2O remains as a gas; in 
a right temperature band, and relative temperature fluctuations, we get a cloud; 



Albedo 




ThermaUzation - 



- Greenhouse 
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hence albedo. Clouds have double effect: reduce the greenhouse because they steal 
H2O molecule from the gas phase, moreover reduce by albedo the incoming radiation 
reaching the ground. Less ground temperature means more super-saturation of H2O 
in the air, and rain, and with it latent heat released. But this is only one step in 
the global water cycle. Water in bulk travels and evaporates elsewhere feeding again 
the greenhouse job. This complex circulation of mass and energy is self-stable. The 
triple feedback is a powerful governor for both the inorganic dynamics and the life 
dynamics. In conclusion the transport of the asymptotic molecules must be gaseous 
and should also be liquid. 

The feedback belongs to the realm of control dynamics and is, by definition, 
opportunistic. Therefore in this paper which is addressed to extrapolations and 
predictions regarding unknown combinations star-planet we cannot take into consid- 
eration opportunistic issues. 

Bearing in mind these recommendations, we can examine, in the next section, the 
behaviour of the temperature of the protagonist planet, which is assumed without 
atmosphere and consequently without pressure. This theoretical planet should be 
used as a theoretical input, rather than a description of the host living planet. 
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5 Temperature oscillations and total angular momen- 
tum 



We assume for the sake of simplicity that the system star-planet is isolated, and 
consequently, that the angular momentum 

J = L + B 

where L is the orbital angular momentum and a is the proper planet spin, is a 
constant vector. The radiation flux impinging on the surface of the planet at a point 
P, depends on the kinematics star-planet; and the constraint J = const, alone is 
compatible with several possible configurations in coordinate space, each depending 
on the initial condition in the star-planet phase space. Nobody gives the initial 
condition for any two body astronomical system. We should rather think of the 
evolution of each planetary system and consider that on the long time scales the 
evolution can set itself in some sort of near stable dynamical configuration. These 
global configurations may be compatible with planetary elliptic orbits. This happens 
with the solar system, and it may be a general feature. 

We consider our protagonist planet in this regime and we add another simplifying 
hypothesis, namely that both L and a are conserved. The particular values of L and 
fj, of course, come with the evolution of the global planetary system. If we do not 
introduce this assumption, there is very little we can say analytically. 

In conclusion, using L and a fixed, we see in appendix B, that the flux 0*" has 
the form 

= i^(V,e)f(0,(/;,V';C,V'o,c^d) (5.1) 



where is the time-angle on the ellipse, < ^/^ < 2 7r (fig. [A.l| ); 
9 is the colatitude on the planet's surface < < vr; (fig. |B.2D 
ip is the longitude < < 2 7r; 

Wrf = ^ is the pulsation of the spin motion with given period (the index d stands 
for day). 

A given meridian rotates about the fixed spin axis with fixed angular velocity uj^ 
(p = iOdt (5.2) 
Notice that t is related to ip by the time-angle equation t = F{'ip), ( |A.6| ), so that 

ij = F-\t) (5.3) 
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The dependence (^^) of the orbital angle i/; on t is nonlinear, and implicit from 
(A. 6). We cannot write il) = ojyt except in the e = orbits. The angle ■0 is given 



mod 2tt and the function t = F(ip) is given mod Ty. Accordingly, we can define 
(appendix A) 

27r 



with 

r,, = 2 7r 



a3 



7 {Ms + mp) 



Using the explicit relationship between ip and t, t = F{ip) , we adopt ip as the inde- 
pendent variable in the calculation of 



Comments on todi The orbital motion is periodic with Ty] this is the motion 
of the center of the planet in the center of mass system star-planet. On the contrary 
the motion of a point P (appendix B) is not periodic in general. More precisely the 
motion of P is: 

periodic with period = nr^, if ^ = n 

periodic with period n2 Ty = rii Td, if ^ = ^ 

non periodic, if — = real number 

f ' UJy 

The three cases are indicated in appendix A (fig. |A.4| ). 

Using as independent variable rather than t the difficulties related to ^ are 
carried into the function f, or, more precisely, into the way in which the point P 
explores the domain d - n > (appendix B, fig. p. 31 ). 

This is the general kinematics we are starting with. Given the input the 
general non equilibrium thermodynamics of the planet follows. In general the fiuid 
circulation on the surface layer of the planet has three main inputs: 

- thermal exchanges with the solid surface which receives the radiation and ther- 

malizes it, 

- direct interaction with the radiation in the form of albedo, 

- heating of the atmosphere due to greenhouse 

Notice that the albedo in its turn is generated by underlying phenomena in the 
thermalization domain, namely phase transitions either permanent (ice) or going on 



16 



and off (clouds). Consequently there is a complicated thermodynamical loop rather 
similar to a control feed back. 

The organic ecosystem participates to this feedback. So both phase transitions 
and global life are in some sort of "dialogue" with the star. This issue is outside 
the scope of this paper. Here we limit ourselves to the non equilibrium thermody- 
namics response to the stellar input in the first approximation, namely the Fourier 
equation (see Appendix C). Notice that if we imagine a perturbative expansion of 
the equations governing the interaction between the fluid component and the rigid 
component, the Fourier equation would be clearly the first order approximation. In 
this approximation, the temperature field follows the input linearly, except for a time 
delay §. 

After this long introduction, we return to the spin and orbit periodicity, namely 
ojd, Notice that if ni is large and n2 is small, the spin typical time is, in a 
certain sense, the dominating clock which misses the period Ty by an error which is 
a little percent of Ty. So the important clock is r^. This happens with </>'", and with 
T as well. 

We know that in the terrestrial realization of life the important clock is r^. 
Animals and human beings act along the time span of a day with accuracy of seconds 
or minutes. On the contrary the phenomena recurrent with the year are met in 
general with inferior accuracy. Clearly living organisms are driven by the daily 
attr actor. 

For a planet having a combination and cOy, with ni and n2 both small, neither 
Td nor Ty can be good clocks. This poses difficult problems and, as terrestrial animals, 
we are unprepared to face them. 

We limit ourselves to present a set of figures where we compare a planet that has 
the same ( ^ 23° as the Earth and the same value of the eccentricity, e = 0.017, but 
with 

Wrf = 6 

instead of the canonical uJd = 365 ujy . 

For this low value of e there is quasi symmetry between the north and the south 
hemispheres. In other words north and south are simply a way of labelling the 
stars that a terrestrial observer can see. For large e north and south are not at all 
symmetric, therefore their meaning is no longer arbitrary; they are in fact thermo- 
dynamically different. 

Furthermore low means small polar cups. Remember that the definition of the 
polar cups is the following: 
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0<^<C North ^^-<^<^^- 
0™ = 0; IpAE < < ^SE 



. . e = 0; ^5i? < ^Z' < ^/'AB 

vr - C < ^ < vr South 

0™ / 0; IpAE <^p<^pSE 
SE: spring equinox, AE: autumn equinox. 

The region C < ^ < — C contains the oscillation of the incoming flux with pulsation 
ujd- This is the region where the daily clock operates. 

Figures labelled with a), show the input figures with b) show the daily 
average input cp^^ 



1 

2^ 



2tt 

dip (5.4) 







Figures labelled with c) show the temperature field for the rigid black planet, nu- 
merical solution of the Fourier equation (appendix D). For e = 0.017, we see that the 
equinoxes (indicated with E) and solstices (indicated with S) are almost equidistant. 

In these figures the daily oscillations are so dense, that are not distinguishable 
in the graphs. We have indicated with light grey when the input and temperature 
field contain daily oscillations. Dark grey as in fig. a), a), a) means no 
daily oscillations. 

To familiarize with the graphs showing the calculated temperature field of the 



black planet, we present in figures 5.3 the case ( = 23°, e = 0.017, LOd = 365 



namely the Earth. These figures refer to 6* = 0, north pole, 6 = j, and 6* = f , equa- 
tor. The south hemisphere is not shown since for e ~ it is the complement of the 
north hemisphere. 



Fig. 5.1 shows the north pole. Fig. a): in one year the flux is like one day 
(summer, the grey bump) and one night (winter, the flat = interval). The 
length of the day is equal to the length of the night. The temperature is shown in 
fig. c). Evidently, the non zero temperature during the time in which the radiation 
is absent is due by Fourier conduction flowing from the 9 > Q zone. The small 
temperature oscillations during the time Td are also originated from the 9 > C zone. 
Our graph helps understanding the role of the water phase transitions on the thermal 
behaviour of the planet, comparing this graph with the observed temperature. 



Fig. 5.2 shows 9 = j. This latitude is for the Earth the temperate zone. We 
see that both input and temperature have during the year two dips and two bumps. 
They are what by tradition we call four seasons. The reader living at this latitude can 
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look at the graphs and appreciate how much, and in which way, albedo, greenhouse, 
wind and rain affect his own surrounding thermodynamics. 



Fig. 5.3 shows the equator. 



We move then to the comparison with the case ( = 23°, e = 0.017 but ui^ = Gu>. 



instead of uJd = 365 ujy. Fig. 5A shows 6 = j. We see first of all that the concept 
of season is hardly visible, by considering fig. a). The temperature has a huge 
amplitude of oscillation and a characteristic difference between the cooling branch 
and the heating branch connected to a stronger time delay between temperature and 
radiation input in the descending branch. This property has been discussed in |^. 
See also appendix C. Greenhouse and air circulation in general act in the direction 
of depressing the oscillations that are proper to the Fourier equation. 

In general, we can say that both a high daily frequency and the fluid circulation 
act together as a mixer that homogenizes the temperature field on the surface. 

In conclusion we need acrobatic hypotheses on the fluid circulation to modify the 
above oscillations which are incompatible with the request of Tb lying in a reasonable 
band. See fig. [s!^ c). 

Comments on C, e We keep for the semi-major axis a the same value that we 
have for Earth and we keep also uJd = 365 ujy in order to have something to compare 
with, when we study a different value of C and e. 

Notice that same a means same Ty, and same ratio ^ means same mixing power 
of homogenization on the temperature filed, as discussed above. In the following 
figures we illustrate the case C = 67°, e = 0.5, ujd = 365 . First of all with e = 0.5 



the symmetry north-south is completely lost. Solstices and equinoxes are unequally 
separated from each other. There is a long interval between spring equinox and 
autumn equinox and a short interval between autumn equinox and the next spring 
equinox. 



Fig. |5.5| represents north pole: a long warm summer and a short extremely cold 



winter. At the south pole we get instead (fig. 5.6) a very long cold summer and a 
short hot winter. 

Notice that on the Earth the populations living in the South hemisphere are used 
to call "summer" the period between the autumn equinox and the next spring equinox, 
namely the interval of time that the populations living in the North hemisphere call 
"winter". Here we do not adopt this inversion of the names. 

Moreover for (" = 67° the polar cups are large, the temperate zone is very narrow. 



Let us consider the equator, fig. 5.9. The seasons appear in a wild way. We 



have 4 seasons, namely two dips and two peaks with a long cold summer and a short 
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interval in which are packed the two equinoxes and the winter solstice. 

The middle latitude 9 = j and 6 = ^ are shown in fig. 5.7 and 5.8. These 
latitudes that are temperate zones for the Earth; in this case instead are very close 
to the behaviour of the poles. The temperate zone in fact is a narrow band around 
the equator. 

In conclusion the orbital parameters C = 67°, e = 0.5 are highly unfit for life, 
unless acrobatic hypothesis on the atmosphere are injected into the general picture. 
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Figure 5.1: ( = 23°, € = 0.017, ujd = 365 cjy, ^ = 0. North pole. 
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Figure 5.2: C = 23°, 6 = 0.017, LJd = 365 e = ^. 
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Figure 5.3: C = 23°, 6 = 0.017, LJd = 365 Uy, 6 = ^. Equator. 
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Figure 5.4: ( = 23°, € = 0.017, Ud = 6oJy, 9 = f . 
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Figure 5.5: ( = 67°, 6 = 0.5, ujd = 365 cjy, ^ = 0. North pole 
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Figure 5.6: ( = 67°, € = 0.5, LUd = 365 U)y, 9 = TT . South pole. 
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Figure 5.7: ( = 67°, € = 0.5, iOd = 365 iUy, 6 = j. Note: the meaning 
of the darker area in a) is explained in the magnified drawing on the upper part of the figure: the 
incoming radiation oscillates with pulsation uia, but there is not radiationless night 
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Figure 5.8: C = 67°, £ = 0.5, LJd = 365 UJy, 9 = ^. 
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Figure 5.9: ( = 67°, 6 = 0.5, ujd = 365 cjy, ^ = f . Equator. 
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6 Concluding remarks 

We have illustrated in §5 the cases a;^ = 365 tOy and uJd = Qujy- We did so for 
ease of calculation and because the year's and day's cycles indeed dominate the 
surface thermodynamics of the planet in the global stationary regime. Choosing two 
examples satisfying uja = nujy is therefore significant. 

On the other hand we remark that the terrestrial ratio ^ is not exactly 365, 
and in general the rule ^ = in. cannot be expected to be the law of nature. In fact 
complicated astronomical events are responsible for such ratio. 

When ^ = real number, the temperature field is more complicated and we have 
shown in appendix C that the average field < T >e,ip (t) has annual periodicity, in 
the stationary regime, but the field T{R, 9, if, t) is non periodic. This means that 
the Fourier expansion of T(R, 9, ip, t) spans the infinite interval 

< t < oo 

and therefore T contains a continuous spectrum in the frequency ly. Instead the 
average < T >e,ip {t) spans the finite interval 

t' <t <t' + Ty 

and therefore < T >g^^ (t) contains a discrete spectrum: 
_ ^ 

This strange fact comes as a consequence of the simplicity of the elliptic motion, 
which determines for the radiation input the factorized form 

r = K{t)<^{9,ip,t) (6.1) 

Evidently, a planetary motion sensible to many body interactions, produces a radia- 
tion input which is much more complicated than ( |6.lD ||8| and the temperature field 
T{R, 6, ip, t) is extremely difficult to analyze. 

The above comments come from the fact that the equation of motion (see ap- 
pendix C) 

dT 

CvPv-g^ = K \;e,v T + {9, ^, t) 5{r - R) (6.2) 

is linear dissipative and is attracted by the driving term ( |6.lD . 

We can make at this point a step forward and move to the consideration of the 
fluid motions on the surface of the planet, see §4. The partial differential equation 
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of the fluid motions, the Navier-Stokes equation containing the heat transport, is 
non hnear, and is coupled to the linear equation ( |6.2| ) which acts as the free field, 
or unperturbed term, in the self-consistent description of the interaction fluid-rigid 
surface. The fluid motions are driven by the field T solution of (|6.2[ ), and in turn 
operate on it according to some iterative procedure, not discussed in this paper. 

Now, the fluid velocity field has the property of being chaotic. If we consider the 
global solution, we are induced to think that it will manifest a regular part and a 
chaotic part in some sort of proportion. This concept of "proportion" is the extrap- 
olation of the property discussed in this paper: periodic the surface average, non 
periodic locally. It is exactly this proportion which is the most important constraint 
that the planetary orbit induces in the relationship Star-Life. 

The living organism must cope both with the regular, periodic part of the am- 
bient, and with the remaining chaotic unpredictable part of the ambient. Therefore 
its control strategy must be quite complex, perhaps we may say intelligent, and able 
to deal with predictable ambient inputs plus unpredictable ambient inputs. This 
analysis could be rejected in favor of the principle that unpredictability is dictated 
by geography. In this view the discussion of extraterrestrial life seems to be re- 
duced to casualness, so that we return to the geocentric viewpoint. This is not the 
viewpoint adopted in this paper. Above the uncertainties of geography stand the 
general constraints discussed in §5, where we show that there are acceptable or non 
acceptable combinations of star temperatures and orbit semi-axis, plus acceptable 
or non acceptable combination of eccentricity and inclination. The definition of con- 
straint we have adopted rests solely on general thermodynamics, and this is the only 
way we have to talk of the compatibility with life on planets which are not directly 
observable. 
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Appendix A. The motion on the elUpse 



Let us consider the ellipse of fig. A.l. This drawing corresponds to eccentricity 
e = 0.8. In this drawing it is also indicated the projection of the spin a in the orbital 
plane, a a-, ipo is the orbital precession. At point F (the first focus) there is the star. 
At point P there is the planet. 



— — 0--' 




Figure A.l: 



d, ^ are polar coordinates centered in F. 

d(ip) = PF is the distance star-planet. PiF = perihelion; 

P2F = aphehon; 

PiO = P2O = a = major semi- axis; 

FO = F'O = c = coordinate of the focus in the Cartesian system centered in O; 
c 

e = - = eccentricity; 
a 

OQ = b = a Vl — = minor semi- axis; 



RF =p = a{l-e^ 

The ellipse equation in the above polar coordinates is: 

d V = -— = --, ^ (A.l) 

1 — e sm ip 1 — e sm ip 

Prom (A.l) we see explicitly the dependence on the parameters a and e. In gravi- 
tational dynamics the physical parameters are Mg, star mass, nip, planet mass, 7, 
Newton constant. The equation of motion is: 



The solution of the vector equation of motion ( [A.2| ) is a two dimensional trajectory 
lying in the plane orthogonal to L = fid A v (the ecliptic plane). The trajectory 
is determined once we fix £^ = £'p + V(d), total energy, and L =\ L \. The two 
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parameters a and e, and the two constants of the motion E and L, are related by 
the following equation: 



E = -7- 



2a 

(A.3) 



L = J a n {1 - e"^) nip Ms 

Since the motion obeys a time reversal invariant equation, the origin of the time is 
arbitrary. We set V'o = and adopt this convention relating and t 

tp = spring equinox t = SE 

■0 = f summer solstice t = tss{e) SS 

tp = n autumn equinox t = tad^) -A-E (A-4) 

ip = 3^ winter solstice t = t«,s(e) WS 

tp = 2-K spring equinox t = Ty SE 

These names are related to the fact that for ip = the projection of the spin a on 
the ecliptic, ffa, is perpendicular to d: 

If we take into account the precession ijjQ, we get 

(<7a,d) = - + 

Keeping always the convention i/j = 0, t = 0, we have in the case of precession ipo 
^ = ipo spring equinox t = to(V'o) SE 

■i/' = f + V'o summer solstice t = tss{e) SS 

tp = n + ^0 autumn equinox t = tae{e) AE (^-5) 

ip = 3^ + ^0 winter solstice t = tyjsie) WS 

ip = 2n + tpo spring equinox t = to + Ty SE 
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(see also fig. A. 2 ). Notice that '0 is the natural parameter that describes the elliptic 
motion. On the other hand the time t is the natural parameter for the description 
of the physical and biological dynamics. 
The relationship between tp and t is given by 



L 



p\ij)d^ = F{ij) 



The explicit calculation gives 



-2a^ y a/1 - e2 tan"^ 



e — tan ■ 



For e = we have the result 



. 2 ^ , 
t = a -ip; 



and 



L 
a? fi 



2 fi e{l-e'^) cosil^ 

+ a — ■ h const (A. 6) 

L e sm tp — 1 



(A.7) 



The index y stands for year. 



SE 



¥o=0 



SS 



AE WS SE 



"Summer" 



"Winter" 



WS 



"Winter" 



\\fQ- Till 



SE 
\|/=7r/2 



SS AE WS 
\|/=7i: \|/=37i:/2 \|/=27t 



"Summer" 



"Winter" 



Figure A. 2: T/ie relationship time-angle e = 0.5. 



From ( |A.6| ) we get the annual period 

Ty = + 2 vr) - t(^) = 2 TT ^ Vl - = 2 vr , 



L 



7 (Ms + mp 



(A. 



Using the time-angle equation (A. 6) we can calculate the time intervals separating 
the typical values labelled solstices and equinoxes. Notice that, obviously, for the 



34 



circular orbit we have (| 



2 

t = a Y il^; or ^ = ujyt 



Moreover the relationship time-angle is linear, therefore the solstices and equinoxes 



are equidistant in time. In fig. A. 2 we show what happens when e = 0.5, both for 

1^0 = and ^^0 = 1 

The time-angle equation ( |A.(]| ) is important (see appendix B). We show in fig. 
[A.3| the relationship ip, t in the case of e = 0.017 (this is the terrestrial orbit almost 
circular) and the case of e = 0.5. 
A point Q on the surface of the planet with polar coordinates R, 9, ip, has a 




Figure A. 3: 



complicated motion in the reference system of the star. We give a simplified example. 

If the spin of the planet is a constant vector perpendicular to the ecliptic plane, 
the point Q, at fixed 9, rotates around the spin axis with angular velocity uJd {d for 
day). We can write 

ip = Uldt 

The ratio between u)d and ujy = y- can be a rational number or a real number. These 



cases are illustrated in A. 4 
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niujd = l iOy 




Ud = real number • 



Figure A. 4: Motion of a point Q on the planet surface. The point Q rotates around the spin 
axis located in P with angular velocity ujd- The point P moves on the ellipse with focus F, position 
of the star, with angular velocity not constant. The relationship '-^ t is given by In the 

figure we have chosen e = 0.5. 
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Appendix B. The incoming flux 0*"^ 
The Planck diluted spectrum 

Let us discuss first the concept of diluted Planck spectrum. We refer to fig. 



B.l 




The conservation of energy in the solid angle dO, gives us 

RlaT^ = (f af^ (B.l) 

T depends on T^, Rg and d, the distance star-planet. Remember that for elliptic 
orbit d = d{ilj) and T = T{'ip). 

T is the thermalization temperature, on a sphere of radius d, of the photons 
emitted by the star according to the Planck law 

l{T,,v)dv= '^^Y dz.; \ZM = ^ (B.2) 

- 1 ™ 

h is the Planck constant, c is the speed of light and k is the Boltzmann constant. 
The spectral distribution ([B.2|) at a distance d is diluted by a factor Vg. 



T>2 rpi 

In other words at the distance d arrives the diluted spectrum: 

IduiTs,i^)=Vs-IiTs,i^) (B.4) 

The black surface positioned at d has temperature T and in its turn emits with the 
spectrum 

T(f , u) du = ^^-^ du (B.5) 

ekf — 1 
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The dependence of 0'" on tp, 6, Lp; ( 



We can write the incoming photon flux (p™' for a point P on the surface of the planet 
with polar coordinates R, 6, if, where: 
R radius of the planet, 

< 9 < n is evaluated starting from the north pole {6 = 0), and ending at the south 
pole 6 = TT. 
< V9 < 27r 

The point P, at fixed 9, rotates around the spin axis with angular velocity ujd {d for 
day). We can write 

and is important to establish the relationship between = 0, t = 0, ^|J = 0. We 
chose 



ip = 0,t = 0, ip = to be noon at spring equinox 

With this convention we have 



cP'''ii;,9,ip,t)=aT^ 



R„ 



d^iiP) 



d-nQ{d-n) = K{i;)f{9,^,^lj) 



using ([A.lp 



f)2 n2 



(1 — e sin?/;)^ 



(B.6) 



(B.7) 



(B.8) 



(i2(^) " ^ a2 (1 - e2)2 
i\ w 

dS is the flux (in W) impinging on the surface element dS positioned at P and 



notice that dim = ^ 



having orientation n as shown in fig. B.2 



To understand eq. (B^), consider fig. BT and B^ jointly. The density of flux 
is the radiation seen from the planet, where it is the star that rotates according 



to the angle ip, (fig. B.2). The unit vector n is the direction of OP, where OP has 
coordinates R, 9, ip. The unit vector d is directed from the center of the planet O 



to the focus F (center of the star), fig. A.l and B.l 



Notice that in fig. \B.£ the center of the planet is called O, while in fig. A.l the 



center of the planet is called P. 

The unit vector d has direction opposite to the momentum of the incoming pho- 
tons 



d 



1^7 1 
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Figure B.2: 



After these remarks on the notation, by consideration of fig. B.2 we obtain 

d ■ n = cos 9 ■ sin(^ — V'o) " sin Q + sin 9 ■ cos{ip — V'o) " cos{ipo + a;^ t) + 
+ sin 9 ■ sin(^ — ^o) " cos ( ■ sm{ipQ + uj^ t) 



(B.9) 



ipQ is the precession. In the following we put ipQ = 0. ipQ is a given meridian. 



According to the origin of time (|B.6D we put ip^ = 0. There is no special identification 
for a meridian. 

The ecliptic is the plane i' , j' with normal unit vector k'. The obliquity C is the 
angle k, k' . d rotates in the plane i' , j'. The reference system of the planet is i, j, k 
The unit vector n rotates around k. The unit vector d is first rotated from the initial 
position d = i' to the current position d^ip) with a rotation Rk'- The unit vector n is 
carried to the reference system i' , j' , k' with a rotation R(^. From these observations 
follows that the scale-product d - n \s given explicitly by 



\ 



/ 1 

COS C sin Q 
y — sin Q cos C, j 



Ri 



I cosip siml^ ^ 
— sin ^ cos ip 



V 











1 



d ■ n 



10 



^ cos^ sintp ^ 
— sin ip cos ip 
1 



/ 1 
















COS C sin C 
— sin C cos C J 



sm tl cos If 
sin 9 simp 
cos 9 



39 



which is (|B.g| ). 

0(x) is the step function: 



e(x) 



1; x>0 
0; X < 



The expression d • n is a function of three variables ip, 6, (f, therefore the constraint 
d-n> 

determines a three dimensional domain in the space ip, 9, ip. 



We show in fig. B.3 the domain {0, ip) for various values of ip. In these figures the 
precession is set to the value ■00 = 0- The parameter C, affects the domain d ■ n > 
in a very interesting way: we show the dependence on presenting two columns, on 
the left C = 23°, on the right C = 67°. 

Harmonic expansion of 0*" 

Remember the time angle relationship ( |A.6| ) 

can be written as 



= K{i;) f{e, ip,^P)=a Tt ^ (1 - e sin V)' f (^, p, V) (B.IO) 

(1 — e^)^ 



with ip = LOdF{'ip). Alternatively 0*" can be written as 

0- = Kit) f{9, ip, t) = aTt (1 - e sin f{e, p, t) (B.ll) 

(1 — e^j^ 

with tp = F^^{t), p = iOdt- The expression ( |B.1C| ) is more natural considering the 
explicit expression of d • n and the fact that t = F{ip) is given explicitly in ( |A.6| ). 
Vice versa ip = F^^{t) is given implicitly, or numerically. 

We have seen that Kiip) is periodic in the time-angle ip because is periodic. 
The factor f(0,99, 0) requires a special discussion. Consider eq. (|B.g| ). We see that 
Ip enters in terms smip and cosip, which are periodic 2tt, and enters also in the ip 
dependence, which is of the kind sin ip and cos ip. 

Since 99 = u;^ t we write 

ip = L^dF{lp) 

The question is whether cos{^iOd F{ip)) or sm(LOd F{7p)) are periodic in ip or not. The 
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Figure B.3: 



41 



answer is yes if uj^ = nujy, in fact in this case we write 

2 vr 

UJdt = nujy F^ip) = n — F[il)) 



but 



r, = F(2^)-F(0) 



and we can put -F(O) = 0, using the integration constant in ( A.6 ), therefore 



ijjdt = 2 Tin 



F{2tt) 



and consequently cos ( 2 vrn ^f/jtr 1 ) periodic with period riTy. On the contrary 



when ujd = ctiOy, with a real number 
F(7p) 

Lu^t = 2tt a — r 

F{2tt) 

In this case uJdTy is no longer a multiple of 2 7r and in other words f(0,99, ^) is no 
longer periodic in ip. This remark is very important for the study of the field T(f, t) 
which is driven by cj)^'^. 

In the following we discuss the case uJd = nujy. 

Let us consider the harmonic expansion in namely given by (plol ). The 
harmonic expansion of (B.IO) utilizes the orthonormal complete set 

1 



'2^ 



(B.12) 



tM 



(B.13) 



/'2 TT pItt 

dT= dip sin^de* 

8 7r2 Jo Jo Jo 



The expansion is 

£,m,M 

^ = 0,1,2, ... ; m = -£,-i + l 



? ' ' ' ^5 



with 



imM 



I 

Jstt^ 



dTrw^:{e,^,^) 



M = 0,±1,±2, 



(B.14) 
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Notice that are complex, 4)irnM are complex, but (jf^ is real. We know that 



4 vr / 4 vr 4 



so that 



M Lm 



(B.15) 



(B.16) 



but 



- V / doy,^ = yoo = ^ 



so that ( B.16 ) becomes 



M il> lin 



M 



Since K['iIj) is known explicitly (eq. 



K{^)=aT: 



' a?{l- e2) 



2\2 



1 — 2 e sin ^ + e 



namely 

2 1 — COS 2^ 



2 1 — COS 2^ 



P(e) 1 - 2e sinV' + e 



Comparing like terms in ( |B.17| ) and ( [B.18| ), we find 



im 
''00-2 



vr 



4^/2 



im 
'^002 



vr 



4^2 



62P(6); 



(B.17) 



(B.18) 



*^001 



vr 



(B.19) 



"ooo 



vr 



2^/2 



(2 + e2)P(6) 



Obviously the component ooo is the largest and the only one that survives in the 
e ^ limit. The annual average of the incoming flux is 



(B.20) 



Next let us consider the harmonic expansion in t, namely 0*" given by ( B.ll| ). The 
factor K[t) is periodic in t with period Ty. The orthonormal set is 

1 i^Mt 



Yem{0,(p) ■ e 



fB.211 
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The expansion is 



l,m,M 



(B.22) 



0,1,2,...; m = -£,-£+!,...£; M = 0, ±1, ±2, . . . 



with 



J Ty J 

Let us perform the average over 9, 



2^ Mt 
y 



(B.23) 



But we know that 

< >e,^=lK{t) = ^ KMe'^""* (B.24) 



4 ^/t\. 



The expansion ( p.23| ) contains all the Fourier components, M = 0, ±1, ±2. . . and 
is not restricted to five Fourier components as it is with the expansion in ip, formula 
(|Bl9D . 



Concerning the non averaged expansion ( p.22| ) in the case we are considering, 
namely uj^ = nujy ox Ty = n we must have a component M = n that is the second 
large mode of the expansion. Nevertheless eq. ( B.22| ) contains all the other modes. 



It is clear that the expansion in ip is helpful in order to appreciate the kinematics 
of the elliptic orbit; the expansion in t is necessary for the study of the Fourier heat 
equation given in appendix C. 

Finally we consider the case ^ =real number. In this case the domain 9, (p is 
always bounded by 47r, but the domain t is unbounded. It follows that the ortho- 
normal set corresponds to discrete indices with respect to 0, (/? and a continuous 
index with respect to the time t. 

The eigenfunctions are 

iy/^ = -Ly,„(0,<^)-e'-* (B.25) 
y It: 

and the orthonormality relations 

— dn dtY;,^,Ye^e'^''~'''^' = 6,,e'Sm,m'Siiy-u') (B.26) 

^TT Jo 



44 



The harmonic expansion is 



J2 J di^<PZu-wLio,v,ty, 

£,m 

= 0,1,2, ... ; m = -£,-£ + 1,. 



(B.27) 



with 



POO 



The interesting fact is that (p'^'^ is non periodic in t but its 9, ip average is periodic 
in t. 



1 1 



y 



as we have seen above in formula (B.24). 

Dimensional remarks for the expansion in t. 

. W 



dim I 



m' 



(B.28) 



dim{W,'iie,^,t))=s 



dim {Wf'^^{9,ip,t)) = dimensionless 



W 1 

dim{(t)ZM) = — S2 
m^ 



W 

dim{(pZ^) = — - s 



m' 



dim = s 
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Appendix C. The temperature field. 



The study of the temperature field develops along these steps. 

1. The planet is assumed to be homogeneous, namely density, c^, specific heat, 
K, Fourier coefficient are constant in the crust. The planet is made of a hot 
core and a crust according to this figure. 




For r < Rq there is a hot core which cools off and produces the source for the 
endogenous temperature field decreasing from Rq to R. 

From the star comes the external input For r < Ri the temperature flow 
generated by the input is negligible. With the assumption of homogeneity 
the Fourier equation for the temperature field T(r, t) is conveniently written 
in polar coordinates 



dT 

Cv Pv = n^r,e,>fT; Ro<r<R (C.l) 



^ _ 92 2 5 f r ^2 Id 1 92 

dr^ r dr | dO"^ tan 9 89 sin^ 9 dip"^ 

Dimensionally eq. ( |C.fD is As equation ( |C.1| ) is linear, the general solution 
is the sum of the endogenous solution belonging to the boundary condition at 
r = Ro 

T{Ro,t)=T,ore{t) (C.2) 

the cooling of the core, plus the solution driven by the incoming flux (p^'^{t). 
Both and Tcoreit) are inputs independent of each other. We write 

Tgen{r, t) = T,ore{t) + T{r, t) (C.3) 
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The behaviour of the temperature field from Tcore{t) has been studied in a a 
preceding paper Here we concentrate on the star driven temperature field. 

The statement of the problem implies 

- The differential equation ( p.!] ) 

- The imposed boundary condition, or driving term 

- The initial condition 

We can take into account the driving term writing 
dT 

= KAr,e,^T + (l)'''ie,ip,t) 6{r - R) (C.4) 
with 5{x) Dirac delta function. 

2. Let us comment the dissipative nature of equation (|C.4| ). Consider a shell 
r — 5 < r < R. (in the following we will justify the concept of shell 5) and take 
uniform temperature within this shell 

T{R - 6, e, if, t) ~ T{R, 9, if, t) = T{t) 

with initial condition 

T{t = 0) = To (C.5) 

Correspondingly, consider the total flux of The index s stands for surface. 
Such quantity has been discussed in appendix B. 

¥^ = ^K{t)R^ (C.6) 
Now we integrate equation ( p. 4]) in the volume 

V = A^{R-~5f~5 (C.7) 

having external surface S and internal surface Si. We get 

c^PvV'^^fi I VT-nd5 + $i"(i) (C.8) 
di Js 

(we have used the Gauss lemma relating volume integral to surface integral and 
we have neglected the contribution to the surface integral across the internal 
surface Si). Now the term 

K / VT-ndS 



s 
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is the outgoing flux. As the planet is positioned in the vacuum, this term is 

VT-ndcS = $""* = 47ri?2aT^(t) (C.9) 

Putting the pieces together we have 

c^p^V—=7rR^K(t)-4.TTR^aT^(t) (C.IO) 
dt 

The non hnear equation (|C.10D has the periodic forcing term 'KK{t)FI?. Let 
us consider the period average of the forcing term 

1 fi+'Ty 



/t + Ty 
K{t') dt' = Ko = constant 



and consider 

rt+T, 



— / T{t')dt' = T{t) 
Ty Jt 

FinaUy consider the segmented equation B 



c^PyV— ^KottR^ -AnR^af^it) (C.ll) 



This equation has the asymptotic fixed point 



4cj 



(C.12) 



Linearizing around Tas, we write 

f = fas + C (C.13) 

and we get the equation 

c^p^Vi = -WTTR^aflC (C.14) 
which has solution 

e = eoe^* (C.15) 

where 

^^_16.Ji'.Tj.. dim(A)=s-' (C,16) 

Cy Pv V 

is the Liapounov exponent of the stable fixed point Tas- 

In conclusion T{t) is asymptotically oscillating with period Ty. We use the 
preceding analysis as a justification to restrict the study of T(r, 9, ip,t) to the 
asymptotic regime, where there is no longer memory of the initial condition, 
because the field T(r, 0, ip, t) has settled down on the periodic attracting input 
term (f)'^^{6,ip,t); in other words we neglect the transient. 
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3. The harmonic expansion is a way to diagonahze the differential equation ( |C.4| ). 
Such program cannot be carried to the end because we do not have the expUcit 
expression of </'£^jy^ , nevertheless from the harmonic expansion of T we can get 
important informations. 

Let us expand T using the orthonormal complete set 



W, 



EmM 



with 



'y JAn 



(C.17) 



(C.18) 



The harmonic expansion is 

r(r, 6, ip,t) = ^ TimM{r) • WlmM 
e,m,M 

1 



y 



y e,mM 

0,1,2,...; m 



' + Af = 0,±1,±2, ... 

(C.19) 



where 



TemM{r) 



1 



dt / dO r(r, 0,99, t)y,% (0,99) e" 



iM 



4 7r 



'y -Jry 

Inserting ( |CrT9| ) into dCTlD we get 
^2 2d 



trriM ■ 



£(£+1) ic^pyMujy 



'-tmM 



(C.20) 



(r) = 0; 
(C.21) 



We project eq. (|C.2l|) on W^,^,}^, and we integrate in dVt dt. Using the ortho- 
normality (C.18) we obtain 



dr^ 



+ - 



r dr 



(C.22) 



4. The radial behaviour. 



Equation (C.22) is a particular case of the Tricomi equation [10|. We are 



interested in the solution of (C.22) which is regular at r = and behaving as 

(C.23) 



a positive exponential for large r; from |1C] we get: 



^ / Ct; Pv UJy 
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Equation ( C.23 ) contains a relevant physical information. In fact consider 
the mode M = 1 corresponding to the frequency ^ and the mode M = n 
corresponds to the frequency ^ = Let us define "penetration length" for 
the mode M the value 6m such that 



exp 
We have 



(R - ~6m) - R 



I 2k 

Cv Pv M LOy 



Cv Pv J. J 



(C.24) 



(C.25) 



So the various higher Fourier components generated by the input frequencies 
Vy = ^ and fd = with large M, have decreasing penetration. In other 
words, they affect a thin layer near the surface. Eq. ( |C.23D tells us which is 
the radial dependence of mode Mi = 1 and the radial dependence of the mode 
M2 = n. From ( C.25|) we see that 

5mi 



M2 
Ml 



(C.26) 



This square root law governs the penetration of the daily temperature oscilla- 
tion with respect to the annual temperature oscillation. 

Notice that when uj^ = aujy, with a irrational number, the Fourier analysis of 
T with respect to the time t is continuous. No matter how complicated the 
spectral analysis of T is [|^, namely the explicit dependence of Hirni/ on u, the 
square root law ( p. 26 ) remains valid. 

5. The angular behaviour. 

Note that the coefficients H^^m appear as a common factor to the addenda of 
equation ( C.22| ), and remain, at this point, undetermined. To find HimM we 
need to use the boundary condition, namely the equation 

=-aT\R) + (C.27) 

In fact inserting the expansion ( |C.19 ) with the explicit solution ( C.23| ) evalu- 
ated at r = R into ( [02^ ), we get an algebraic equation involving Hi^m and 
'^imM- principle this equation gives the coefficients HiraM- Unfortunately 
there are two difficulties. 

The coefficients are not known analytically, and this happens because f(0, (p,t) 
is a discontinuous function. This function could be rounded off with some 
approximations . 
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The second difficulty is that the algebraic equation relating i/^mM and 4>hnM 
contains the fourth power of a series. In this case is necessary to perform a 
truncation. We have studied both procedures, rounding off and truncation, in 
a preceding paper, concerning a particular case where the function f{9, ip, t) 
was simpler. Here we skip the analysis and we leave the study of the angular 
dependence of the field T to the numerical calculation, appendix D. 

6. Finally we consider the case ^ =real number. In this case we adopt the same 
harmonic expansion that has been used in appendix B for the input function 



with 



dt J dnTWl*^i9,ip,t) 



(C.28) 



(C.29) 



The expansion (C.28) is unrestricted with respect to the time interval. As 
we are dealing with a dissipation problem (the Fourier heat equation) the time 
runs from t = to t ^ oo rather than — oo < t < oo as we have for Hamiltonian 
problems. 

On the other hand the solution T of the equation ( p.l|) , as we have seen in 
( C.12| ) is asymptotically attracted by the average driving term in the sense 
that 

1 

1 

a 



<T>i 



(C.30) 



This result tells us that in the asymptotic regime ( C.30I ) we can write 



— / dnr 

Ait 



vr T, 



du Too^ e 

Mt 



iut 



(C.31) 



y M 



Equation ( C.31 ) contains the property that the solution T of equation (|C.l] ) is 
attracted to the periodic behaviour ( C.31| ). 

Another way to express the same result is to write the average of T in the 



interval of time r„ 



t+T, 



V 11 

dil < T >e,^= - 

Ty ^/8^ 

1 1 1 



t + Ty 



dt' / duTooiyC 



iut' 



t— >oo 



t^oo 



Mt' 



M 



_1 ^300 



(C.32) 
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This is in agreement with the fact that the solution of a differential system 
lying on the attr actor has a smaller number of dimensions; in this formulation 
this statement means that the continuous variable u tends to the discrete index 
M. 

Dimensional remarks for the expansion in t. 

dim (T) = K 

dim(l^,^(0,(^,t))=s-5 

dim {W^^„X6, ip,t)) = dimensionless 

dim (T£mM) = K • s2 

dim (T^rnu) = K • s 

dim {5{i')) = s 



52 



Appendix D. The numerical calculation of the field T 



The equation of motion is 
dT 

c,p,-^ = K Ar,e,^ T + t) 5{r - R) (D.l) 

with 

-KVT-n = aT^; r = R 

-KVr-n = 0; r = Ri 

Since the input term is written exphcitly in polar coordinates, this imposes for the 
discrete spatial grid the polar choice. T is then the unknown matrix with three 
indexes, Tm,ij and these indexes are defined as follows: 

<m< mmax] Ri <r < R; hr = ; 

'^max 

2 vr 

0<i<imax; 0<v9<27r; hi = ; (D.2) 



vr 

0<j<jmax; < 6* < vr; h 



] — • ' 

Jmax 

The input term in ( p.l| ) is a given matrix with two indexes: 

namely m is fixed to the surface value m = rrimax i,j run according to (D.2) 



The operators gradient and Laplacian are written in a standard way as relationships 



between respectively first and second neighbours [12|. 



dT 

The step of the time enters in the derivative — and 0*". The time appears 
in (/>™ in two places: 

sin iOdt or cos (J^t 
simp{t) or cosTp{t) 

The function V'(^i) niust be evaluated numerically at each step tj. 
The calculation begins with an initial value 

Tm,i,j{t = 0) = To 

and proceeds with various choices of At, TTimaxy ^max^ jmax- The integration in t is 
calculated with the Runge Kutta method (4th order) The attractor regime is 
found when the value Tm,i,j{v) evaluated after a number of time steps = is 
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sufficiently close to T^_ij(0). 

In practice, for our model, we have used 

f^max ~ 20; imax ~ 30; jmax ~ 40; 

At = 3 minutes; 

c„ = 840J kg-^K"^; p = 2000kg m"^; K=1.5Wm-^K-^ 

For To = 200 K the transient time required to reach the attractor is ~ 5tj/. For 
R — Ri we have chosen the value 

i? — i?i ~ 6 m 

This value is approximatly twice the penetration lenght 6 of the mode M = 1 
corresponding to one year discussed in appendix C. This is the general plan of the 
numerical study. Of course for the study of the daily cycle the penetration lenght is 
smaller and the steps in r are accordingly concetrated near r < R. 
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